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The  excitation  of  the  eigenfrequencies  of  finite  radar  or  sonar  targets,  of  inhomogeneities  in  elastic  materials,  of  geological 
strata  or  of  the  entire  earth  by  the  impact  of  propagating  waves  (of  electromagnetic  or  acoustic  nature,  or  of  ultrasonic, 
elastic,  or  seismic  character,  respectively!  manifests  itself  in  the  appearance  of  poles  in  the  resulting  wave  amplitudes,  as 
described  by  the  Resonance  Scattering  Theory  (RSTl.  In  the  complex  frequency  plane,  these  poles  relate  to  the  ringing  of 
the  scattering  resonances.  In  the  complex  mode  number  plane,  corresponding  poles  are  connected  with  circumferential  or 
creeping  waves.  An  analytic  relation  between  these  two  descriptions  is  indicated  here,  and  a  number  of  examples  from  the 
above-mentioned  fields  will  be  discussed.  We  introduce  the  concepts  of ' Acoustic  Spectroscopy'  and  of  • Radar  Spectroscopy' , 
respectively,  by  exhibiting  the  target's  resonance  frequency  spectrum  in  a  form  familiar  from  atomic  spectroscopy,  in  order 
to  study  the  shifting  and  splitting  of  resonances  'levels'  under  changes  of  target  shape,  and  to  provide  us  with  possible 
solutions  for  the  'inverse  problem'  (i.e.,  determination  of  target  properties  from  echo  properties  -  here,  resonant  echoes). 


1.  Introduction 

Natural  mechanic  or  electromagnetic  oscillations  of  material  objects  take  place  at  well-defined  eigen- 
frequencies,  determined  by  the  constructive  interference  of  internal  and/or  external  surface  waves  of 
standing-wave  type  around  the  objects.  The  scattering  of  external  scalar  (acoustic)  or  vector  (electromag¬ 
netic  or  elastic)  waves  by  such  objects  at  frequencies  close  to  the  natural  eigenfrequencies  exhibits  a 
resonance  character  due  to  the  excitation  of  the  corresponding  eigenvibrations  [1].  In  a  modal  description 
of  scattering,  the  amplitude  is  governed  by  two  variables,  viz.,  the  frequency  variable  and  the  mode 
number  variable.  A  plot  of  the  response  surface’  [2],  i.e.  of  the  scattering  amplitude  modulus  in  three 
dimensions  plotted  against  these  two  variables,  exhibits  a  series  of  parallel  ‘ripples’  inclined  to  either 
axis.  When  intersecting  this  surface  at  a  constant  mode  number,  the  ripples  appear  as  frequency  resonances, 
and  intersecting  at  constant  frequency,  the  same  ripples  lead  to  mode  number  resonances. 
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The  mode  number  resonances  are  due  to  complex  poles  of  the  scattering  amplitude  known  as  Watson' 
[3]  (or  in  nuclear  physics,  Regge’  [4])  poles  whose  residues  represent  the  mentioned  surface  waves 
(creeping  waves),  as  previously  studied  by  Franz  [5].  The  frequency  resonances  (known  in  nuclear  physics 
as  'Breit-Wigner  resonances'  [6])  stem  from  poles  in  the  complex  frequency  plane,  originally  studied  on 
the  basis  of  the  'Singularity  Expansion  Method’  (SEM)  by  Baum  [7],  In  the  latter  context,  they  lead  in 
the  transient-held  case  to  a  Prony  series-type  superposition  of  decaying  sinusoids  which  may  be  utilized 
to  experimentally  obtain  the  SEM  pole  positions  for  target  characterization  purposes.  We  have  shown 
that  the  Prony  series  represents  a  mathematical  concept  which,  when  summed  over  appropriate  subsets 
of  complex-frequency  poles,  synthesizes  individual  surface  waves  which  constitute  the  physical  entities 
[8-10].  Examples  will  be  shown  of  radar  and  acoustic,  as  well  as  elastic-wave  and  geophysical  resonances, 
and  possible  uses  of  the  resonance  concept  for  the  solution  of  the  inverse  scattering  problem  will  be 
pointed  out. 

The  inverse  scattering  concept  refers  to  the  determination  of  properties  of  the  target  from  the 
characteristics  of  the  scattering  echo  response  [11-13].  It  was  shown  previously  that  for  targets  of  known 
(simple)  geometry,  target  composition  may  be  inferred  from  the  location  and  the  widths  of  resonances 
whose  effects  are  carried  by  the  scattered  echo  [14-21].  On  the  other  hand,  target  shape ,  for  a  known 
target  composition,  may  be  inferred  from  the  spectral  scheme  of  the  target  resonances  also,  as  shown 
by  introducing  the  concepts  of  'Acoustic  Spectroscopy'  [22],  or  of  Radar  Spectroscopy'  [23],  respectively. 
Here,  we  plot  the  eigenfrequency  spectra  of  acoustic  or  of  radar  targets  in  a  way  familiar  from  atomic 
spectroscopy,  and  study  the  shifting  and  splitting  of  the  resonance  'levels'  under  changes  of  target  shape, 
thus  assessing  the  sensitivity  of  the  mechanical  or  electromagnetic  eigenfrequency  spectra  of  targets  to 
their  shapes,  which  offers  the  possibility  of  determining  target  shape  from  a  spectral  analysis  of  the 
resonant  echoes. 
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Reference  [24]  gives  expressions  for  the  modal  impedances  Z„  and  admittances  Y„  of  variously  layered 
spheres;  for  a  homogeneous  dielectric  sphere  of  permeability  n  =  1,  one  has,  e.g. 


Z„  =  ~i 


.*(>  X  i/,,  (X  |  ) 


Y„  =  -  i 


k  [Xljn(Xl)Y 

.k_  X  i/n  (X  1 ) 

’fcu  [xi/„U,)]” 


where  *i  =  k  =w/c(,  Ci  being  the  speed  of  light  in  the  dielectric.  For  a  perfectly  conducting  sphere, 
one  has  Z„  =0,  Yn  =  oo. 

Complex  eigenfrequencies  of  the  spherical  target  are  obtained  as  the  roots  of  the  denominators  of  (lc) 
(‘TE  modes')  or  (Id)  (“TM  modes').  However,  for  penetrable  targets  such  as  dielectric  spheres,  or 
conducting  spheres  with  dielectric  coating,  the  resonances  in  the  scattered  echoes  as  functions  of  frequency 
are  sharp  and  narrow,  as  can  be  seen  in  Fig.  1  where  the  radar  cross  section 

<7  =  ^  £  (-l)"(n  +  !)(<!„ -M  (3) 

of  a  dielectrically-coated  conducting  sphere  (outer  radius  a)  is  plotted  vs.  k0a  for  kua  =  37-41.  Accord¬ 
ingly,  the  imaginary  parts  of  the  denominators  are  small,  and  the  real  resonance  frequencies  xY, 


;  »  gft  38  $8  *8  s8  Sja  s? 


a'X  =  5.9  6.0  6.1  6.2  6.3 

=  37.07  37  70  39  33  38  96  39  58 


Fig.  1.  Radar  cross  section  for  a  conducting  sphere  coated  with  a  dielectric  of  dielectric  constant  f  =  6,  relative  thickness  S  =  0.05, 

and  an  outer  radius  a,  plotted  vs.  x  =  kua  =  2ira/A  (from  [25]). 
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are  found  as  the  solutions  of 

1  /  (LZn )  =  Re{[x/t  ^ '  (x  )]/[xh '  (x )]},  TE  modes,  (4a) 

iV'n  =  Re{[^:"U)]/[x/l:"(x)]'},  TM  modes,  (4b) 

with  multiple  roots  labeled  by  /  =  1,  2,  3  . . .  . 

The  expressions  of  (2)  for  the  dielectric  sphere  are  meromorphic  functions  (which  in  nuclear  physics 
are  known  as  Wigner’s  R  matrix  or  ‘derivative  matrix'  [6])  that  may  be  represented  by  Mittag-Leffler 
series,  e.g.  for  TM  modes  [25]: 

l/(iVJ  =  T  I  >  (5a) 

k  i  =  —co  x  x hi 


the  (real)  pole  positions  being  given  by  the  expressions 
Xn,  =  Z„,(Ci/c) 

where  znt  are  the  zeros  of  jn(z),  their  multiplicity  being  labelled  by  I. 
We  may  introduce  the  ‘S-function’  S„M  =  1  +  2 b„  and  find 

i  V  —  -r {2) 

C™  _  r<0»TH  1  1  n  Zn 
•>"  — *>n  ■  y  _  -<1>  ’ 


(5b) 


(6a) 


where 

L <*> / „  ) 

*!■"  =f  ■^r,sa-H(x)^iirn(x)  (»  =  1,2), 

[Xrl  n  (X  /  J 

the  quantity 

r(0)TM  _  [xAl,  2iC_ 

[xh'nu(x)}’ 


(6b) 


(6c) 


being  the  S-function  of  a  conducting  sphere.  Accordingly,  the  S-function 
C™  _  o(0)TM  f  ~  +  itTn  )(i  Y„  ) 

<  ,  .  v/.  w,  v-l» 

\-(<Tn  -lwn)(lYn) 


(6d) 


with  (5a)  being  used  for  (iVn)~l,  exhibits  the  singularity  structure  of  the  scattering  amplitude.  The 
S-function  thus  has  two  kinds  of  TM  poles,  those  given  by  the  zeros  of  [x/tl,"(x)]'  corresponding  to  the 
(broad)  resonances  of  electromagnetic  waves  propagating  over  the  exterior  of  the  sphere,  and  those 
given  by  the  roots  of  (4b)  corresponding  to  internal-wave  resonances,  only  these  being  sharp  and  narrow. 

For  the  latter  case,  one  may  consider  one  resonance  at  a  time  in  (6d),  and  thus  obtain  the  ‘one-level 
approximation'  [25,  26] 


5™ 


TM  >•  /'i  t/’TM 

«(0)7M  X  nl _ nl 

-L^n  TM  .  / :  /  ^  v  nTM  * 

i  x-xn(  +(i./2  )rni 


(7a) 


the  so-called  Breit-Wigner  form  of  the  S-function  [6],  in  which  the  resonance  frequencies  are 


TM  •TM  / 1  v  /  •  TH 

X  nl  Xn,  {K0/k)<rn{xni 


(7b) 


and  the  resonance  widths  are 


,-w 
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Returning  to  the  Mie  coefficients,  one  has,  e.g., 

.  jly  <>(0)TM  ,  TM  ,  -(1/2)/"™ 

bn-i\lsn  U-  )— -ry  +Tl/2)F™ 


sin  f„j 


where  the  first  term  represents  a  series  of  internal  resonances,  superimposed  upon  and  interfering  with 
the  second  term  which  consists  of  a  relatively  smooth  background  containing  only  quite  broad  resonances 
corresponding  to  those  of  a  conducting  sphere,  together  with  a  specularly  reflected  amplitude  [5]. 

The  resonance  frequencies  of  (7a)  or  (7d)  are  thus  given  by  the  (complex)  pole  positions 

xnl  =  x„-(\/2  )r„i  (8) 

representing  the  ‘SEM  poles’  [7]  in  the  complex  frequency  plane.  For  (physical)  real  frequencies  *, 
however,  the  expressions  of  (7a),  (7d)  will  pass  over  resonance  peaks  of  finite  height  as  x  is  varied,  which 
appear  as  the  foothills  of  distant,  very  (in  fact,  infinitely)  high  mountain  peaks  located  at  the  poles 
themselves  [27].  For  the  real-frequency  case,  one  may  instead  look  for  poles  in  the  (complex)  mode 
number  variable  n,  by  expanding  [28] 

Xnt=Xn,t-(n  -nt)x'n,i+  •  ■  ■  ,  (9a) 

and  with  n,  so  chosen  that  x„,i  =x.  This  brings  the  denominators  of  (7a),  (7d)  into  the  form  n  -nt,  where 
n,  =  n,  +  (i/2  )f  ni,  (9b) 

with 

P =  r  n,l/ x'n,t.  (9c) 

Accordingly, 

(9d) 

n-n, 


and  ni  of  (9b)  gives  the  position  of  ‘Watson  poles'  (or  ‘Regge  poles’)  in  the  complex  n  (mode  number)  plane. 

The  above-described  mathematical  steps  leading  from  (7a)  or  (7d)  to  (9d)  may  be  represented 
geometrically  as  follows.  We  show  the  modulus  of  the  resonance  part  of  (7d)  in  Fig.  2  (the  figure  actually 
depicts  the  analogous  case  of  the  amplitude  of  compressional  elastic  waves  scattered  from  a  fluid-filled 
sphere  [29])  as  a  three  dimensional  plot,  representing  the  ‘response  surface’  above  a  plane  with  one 
axis  being  the  frequency,  the  other  axis  the  mode  number  n.  The  response  surface  shows  a  series  of 
parallel  ripples  inclined  with  either  axis.  If  these  are  intersected  at  a  constant  value  of  n,  there  result  the 
SEM  resonances  vs.  frequency;  but  if  the  ripples  are  intersected  at  constant  frequency,  resonances  in  n 
are  obtained  which  are  a  result  of  the  Regge  poles. 

In  Fig.  3,  we  plot  the  modulus  of  the  resonance  portion  of  b„,  (Id),  versus  frequency  for  n  =  1,4,  7 
and  10  for  the  case  of  a  dielectrically  coated  conducting  sphere. 

The  corresponding  resonance  curves  are  the  first  mentioned  intersections  representing  the  SEM-pole 
resonances;  they  are  seen  to  contain  families  of  resonances  whose  members  (termed  ‘Regge  recurrences') 
shift  to  higher  frequencies  with  increasing  n  due  to  the  inclination  of  the  mentioned  ripples  with  either 
axis.  Each  family  corresponds  to  a  surface  wave,  since  asymptotically, 

Pl(cos  0)occos[(j>  +f)0  +4ir]  (10a) 

so  that  (1)  contain  residue  factors  from  the  Regge  poles  which  have  the  form 

exp  i{±(«i  +  2)0  ±4tr-wf}.  (10b) 


Vv\-V'  -  ■  .‘  V'V-7  V  \-.V  ' v. 


*  '  V*  »  *  • 
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Fig.  2.  Response  surface  of  the  resonance  amplitude  for  compressional-wave  scattering  (with  mode  conversion  into  shear  waves) 

from  a  water-filled  spherical  cavity  in  aluminum  (from  [29]). 


This  represents  a  circumferential  (surface)  wave  with  a  linear  propagation  constant  Re(«,  +  \)/a,  i.e.  a 
phase  velocity 

Ci(x)  =  xc/(rti  +  5)  (11a) 

for  the  /th  surface  wave,  a  wavelength 

A,(x)  =  2ira/(«i+l)  (lib) 

such  that  resonance  (n/  =  n),  n  +  \  wavelengths  fit  the  circumference,  and  an  attenuation  exp(-(9/0()  where 

6,  =  2//\,,  (11c) 

which  is  due  to  radiation.  Figure  4  shows  ‘dispersion  curves’  ci(x )  vs.  x  for  the  TM-mode  surface  waves 
/  =  1-4  on  the  coated  sphere  of  Fig.  3  (solid  curves),  and  on  a  dielectric  sphere  (e  =  6)  of  equal  radius 
(dashed  curves).  The  above  analysis  has  thus  led  to  a  conceptual,  as  well  as  analytic,  connection  between 
SEM  poles  on  one  hand,  and  surface  waves  (represented  by  Regge  poles)  on  the  other. 


3.  SEM  poles  and  surface  waves  in  acoustic  scattering 

The  poles  in  acoustic  scattering  amplitudes  can  be  treated  in  an  entirely  analogous  fashion.  The  total 
(incident  and  scattered)  pressure  amplitude  in  a  fluid  medium  ambient  to  a  solid  sphere  may  be  written 
as  [1] 

P  =  e-l“'  I  \n(2n  +  \){jn(kr)  +  (Bn/Dn)h'"(kr)}Pn(cos6),  (12) 

n  -0 


\V 


J  J 
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where  B„  and  Dn  may  be  explicitly  given  as  3  x  3  determinants,  containing  spherical  Hankel  and  Bessel 
functions  whose  arguments  are  x=ka(k  =  cj/c ,  c  being  the  ambient  sound  velocity),  and  Bessel  functions 
with  arguments  kda  and  ksa  where  kd,  ks  are  the  propagation  constants  of  bulk  dilatational  and  shear 
waves,  respectively,  in  the  sphere  material.  The  SEM  poles  are  the  solutions  of 

D„{x)  =  0;  (13a) 

they  are  the  complex  eigenfrequencies  xni  which  we  label  by  /  =  1,  2,  3  ... .  The  Regge  poles  are  the 
solutions  of 

D„(x)  =  0  -13b) 

with  real  x,  i.e.,  they  are  the  complex  pole  values  viix)  which  as  in  Section  2,  represent  the  propagation 
constant  of  the  /th  surface  wave. 

We  have  solved  (13a)  numerically  for  a  tungsten-carbide  sphere  (characterized  by  a  density  ps  = 
13.80  g/cm3  and  by  bulk  wave  speeds  cd  =  6.860  x  10s  cm/s  and  cs  =  4.185  x  10s  cm/s)  immersed  in  water 
(p  =  l  g/cm3,  c  =  1.4760x  105cm/s).  The  SEM  poles  xnl  in  the  complex  X -plane  are  shown  in  Fig.  5; 
pole  families  corresponding  to  a  given  surface  wave  are  connected  by  solid  lines.  Bollig  and  Langenberg 
[27]  give  analogous  results  for  an  aluminum  sphere. 

We  distinguish  here  the  pole  series  corresponding  to  the  first  (/  =  1)  surface  wave,  which  has  an  entirely 
different  appearance  from  those  of  the  other  surface  waves;  in  particular,  it  has  much  larger  imaginary 
parts  so  that  this  surface  wave  is  more  highly  attenuated.  This  wave  has  been  termed  the  ‘Rayleigh  wave' 
[30],  since  in  the  limit  of  infinite  cylinder  radius,  it  goes  over  into  the  well-known  Rayleigh  wave  on  the 
flat  boundary  of  a  solid  half-space.  The  other  surface  waves  (/  3=  2),  are  referred  to  as  the  ‘Whispering- 
Gallery  Waves’. 

Using  (1  la)  at  the  resonance  frequencies.  Fig.  5  can  be  made  to  yield  the  dispersion  curves  of  Rayleigh 
and  Whispering-Gallery  waves.  A  few  of  these  are  shown  in  Fig.  6,  where  the  Rayleigh-wave  phase 
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velocity  is  seen  to  start  from  zero  at  low  frequency,  while  the  Whispering-Gallery  waves  only  exist  above 
an  individual  cutoff  frequency  at  which  their  phase  velocity  decreases  down  from  infinity. 

The  SEM  poles  shown  in  Fig.  5  correspond  to  •internal'  surface  waves  that  propagate  with  speeds 
comparable  to  the  bulk  velocities  of  the  target  material.  In  addition,  there  also  exist  external,  or  F ranz-type 
[5]  creeping’  waves  that  propagate  in  the  ambient  medium  around  the  scatterer,  with  speeds  comparable 
to  the  corresponding  sound  speed.  Their  poles  cannot  be  shown  in  Fig.  5  because  their  imaginary  parts 
are  of  order  unity  or  larger;  hence  they  are  much  more  highly  damped  than  the  internal  surface  waves. 

Apart  from  spheres,  acoustic  SEM  poles  can  also  be  discussed  easily  for  infinite-cylindrical  targets 
[31].  If  a  plane  wave  is  incident  on  such  a  cylinder,  its  wave  vector  making  an  angle  y  with  the  z-axis 
(i.e.,  the  cylinder  axis),  so  that  the  components  of  the  propagation  vector  k=aj/c  are  kx  =  k  sin -y, 
ky  =  k  cos  y,  one  then  has  a  scattered  solution 

pt  =  eik,t~'wM  I  (2-bn0)inTnH»\kxr)cosn6,  (14a) 

n  =0 

with  an  S-  function  Sn  =  2 Tn  4- 1  given  by 

Sn—S, i-r-rj  r  7T777  Ti-T*  (14b) 

Ln  [z  n 


where 


*'"(*) -TP iyr  ('  =  1.2),  (140 

n  n  ) 

and  where 

HT(kxa) 

5 ■  ■  -wr ,14d) 

is  the  S-function  of  a  rigid  impenetrable  cylinder.  For  the  case  of  a  fluid  cylinder  of  density  p'  and  wave 
number  k',  one  has,  e.g., 

Ln  =  P_KoJ'n(k\a)  (He) 

p  J„(kxa) 

(p  being  the  density  of  the  ambient  fluid).  Sn  and  Tn  may  be  represented  in  resonance  form  as  in  (7a), 
(7d),  and  the  internal  surface  waves  derived  accordingly.  They  are  found  to  propagate  over  the  cylinder 
surface  along  helical  paths  with  pitch  angles  (with  the  z-axis) 

y,  =  tan~Un,/akz),  (15a) 

and  with  phase  velocities 


[cos2  y  -Hm/ka)2}1 


(15b)  being  the  cylinder-equivalent  of  the  spherical  surface  wave  speed  of  (11a). 

For  impenentrable  (e.g.,  rigid)  cylinders  ( p'»p ),  the  resonances  in  (14b),  and  hence  the  internal 
resonances,  are  absent  and  only  the  resonances  due  to  external  creeping  waves  exist.  From  (14d),  these 
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are  determined  by  the  roots  of 


H"'\kxa)  =  0, 


(16a) 


which  for  real  integer  n  represent  the  external  SEM  poles  in  the  kxa  plane.  However,  (16a)  has  been 
solved  by  Franz  [5]  for  real  kxa  in  terms  of  the  Regge  poles  in  the  complex  n  -plane,  with  the  result 

nt(x)  =  x  +  (x/6)U3  c'1,/3q,+  ■  ■  •  (x~kxa),  (16b) 


where  q,  =  1.469354  etc.  This  gives  directly  the  propagation  constant  of  the  creeping  waves,  whose 
dispersion  curves  are  thus  found  as 


(16c) 


for  the  phase  velocities,  and 


for  the  group  velocities. 

As  discussed  in  [31],  very  interesting  refraction  effects  occur  upon  the  generation  of  surface  waves  by 
the  incident  wave  since  yt  #  y.  They  consist  in  a  change  of  propagation  direction,  either  in  the  axial 
direction  for  Franz  waves  and  for  internal  surface  waves  with  Ci  «c,  or  in  the  radial  direction  for  internal 
surface  waves  with  c(  s*  c. 


4.  Surface  waves  on  the  earth 


The  example  of  seismic  waves  represents  a  case  of  surface  waves  that  are  not  generated  by  the  scattering 
of  plane  waves,  but  by  the  action  of  point  sources  (e.g.,  depth  charges,  or  earthquakes).  Apart  from  this, 
the  theory  is  exactly  the  same  as  that  described  in  Section  2.  In  particular,  the  phase  and  group  velocities 
of  the  seismic  surface  waves  are  given  by 


and 


(17a) 


cfPU)  = 


c 

drti/dx 


(17b) 


respectively,  where  x  =  ota/c,  a  being  the  earth  radius,  to  the  circular  frequency  of  the  waves,  and  c  a 
convenient  reference  speed  which  actually  drops  out  in  (17).  These  equations  may  be  evaluated  at  the 
known  eigenfrequencies  of  the  earth  xnl,  at  which  one  has  n(  -» n  (integer). 

The  earth  may  undergo  two  types  of  eigenvibrations,  namely,  spheroidal  and  toroidal  (torsional)  ones. 
Surface  displacement  components  (u,v,w)  in  the  ( r,d,<t> )  directions  for  a  spherical,  radially  stratified 
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earth  may  be  represented  as  [32] 


u=  I  Un(r,<o)Y?(6,<fi) 

n,m 

(18a) 

«=ll 

n.m  ' 

,  an(r,cj)dY™} 

tf'nlr.W - + - f, 

1  86  sin  6  8<f>  \ 

(18b) 

n.m 

I  Vn(r,<o)8Y:  ,3^} 

.  .  ,,  Iln(r,ca)  , 

1  sin  6  8<t>  86  J 

(18c) 

and  U„ ,  'f',,  represent  spheroidal  modes  of  vibration  which  are  uncoupled  from  fl„,  the  toroidal  vibration 
modes. 

Gilbert  and  Dziewonski  [33]  have  empirically  determined  a  total  of  1064  eigenfrequencies  x„i  of  the 
earth,  both  spheroidal  and  toroidal,  essentially  by  Fourier-analyzing  observed  earthquake  shocks.  Using 
these  data,  we  obtain  e.g.  the  phase  velocities  of  spheroidal  (Fig.  7)  and  toroidal  (Fig.  8)  seismic  surface 
waves  [34],  plotted  vs.  the  period  r  =  2-ir/tj,  as  solid  lines  connecting  the  data  (points).  These  dispersion 
curves  are  remarkably  similar  to  those  of  Fig.  4  or  6  except  for  now  being  plotted  against  the  inverse 
frequency  variable.  The  curves  are  labelled  by  the  order  /  =  0, 1,  2. . . .  of  the  surface  waves,  and  are 
compared  with  the  theoretically  predicted  phase  velocities  (broken  curves  or  open  circles)  of  the  so-called 
Gutenberg-Bullen  A'  model  [35]  which  consists  of  a  solid  mantle  overlying  a  fluid  core  at  a  depth  of 
2900  km,  with  an  inner  core  at  5200  km. 


Awl 


M 


Fig.  7.  Dispersion  curves  for  the  phase  velocities  of  spheroidal  surface  waves  /  =  0— 10  (as  labeled),  obtained  from  the  empirical 
eigenfrequencies  of  the  earth  (solid  curves  and  data  points),  and  of  surface  waves  /  =  0-9  from  the  Gutenberg-Bullen  A1  earth 

model  (dashed  curves),  plotted  vs.  the  period  (from  [34]). 


In  Fig.  7,  one  observes  a  transition  region  at  c(~8  km/s  at  which  successive  surface  waves  form  the 
dispersion  curve  of  a  ‘pseudomode’,  attributed  [35]  to  the  existence  of  a  Stoneley  wave  on  the  core-mantle 
boundary  which  affects  only  spheroidal  modes.  Our  analysis  obtains  a  second  transition  region  at 
Ci  —  15  km/s  which  is  not  predicted  by  the  Gutenberg-Bullen  model. 
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Fig.  8.  Similar  to  Fig.  7,  for  toroidal  surface  waves  (from  [34]). 

5.  SEM  poles  and  creeping  waves  (transients) 

The  effect  of  SEM  poles  on  transients,  as  well  as  the  case  of  transient  creeping  waves,  will  now  be 
discussed  using  the  example  of  a  perfectly  conducting  sphere,  and  restricting  ourselves  to  TM  modes. 
The  corresponding  TM  part  of  the  first  polarization  function,  (lb),  is  thus 

•  S  ,  2n  +  l  [*/„(*)]'  dP i (cos  d ) 


1  ""  n(n  +  l)[jt/tLuU)]'  dd  ' 

The  complex-frequency  SEM  poles  are  given  by  the  complex  zeros  of  [x/il/'u)]'  which  we  denote  again 
by  xnl.  These  zeros  were  obtained  by  us  numerically  and,  in  Fig.  9  are  plotted  in  the  plane  of  the  variable 
sa/c  =  -ix,  rather  than  of  x  (in  which  the  figure  would  appear  rotated  counter  clockwise  by  90°),  since 
it  is  customary  in  SEM  theory  to  represent  the  complex  poles  in  the  plane  of  the  Laplace  variable  s.  The 
dashed  lines  link  the  poles  according  to  mode  number  n.  When  the  poles  are  grouped  in  (non-symmetric) 
‘layers',  however  (solid  lines  in  Fig.  9,  labeled  by  /),  then  each  layer  will  synthesize  one  (the  /th)  creeping 
wave  [8].  The  present,  correct  system  of  layering  [9]  differs  from  that  of  [8],  and  can  be  rigorously  arrived 
at  as  follows.  The  spherical  Watson  (Regge)  zeros  of  [xh 1  (x )]'  in  the  n-variable  are  given  asymptotically 
by  [5]: 

ni(x)+2—x  +(x/6)l/3  e'^^qi  +  •  •  ■ ,  (20a) 

with  qt  the  same  as  quoted  after  Eq.  (16b).  The  label  /  determines  here  the  /th  creeping  wave,  «/  +  ! 
being  its  propagation  constant,  c.f.  (10b). 

Equation  (20a)  may  be  inverted,  by  replacing  n(-*n  and  solving  for  x  =x„i,  in  order  to  obtain  an 
asymptotic  expression  for  the  SEM  pole  positions  x„(: 

Xni  =«  +2  -(«  +i)1/,6M/3  el,T/\fi  +  •  •  •  .  (20b) 

When  this  expression  is  plotted  in  the  diagram  of  Fig.  9,  the  asymptotic  pole  positions  of  (20b)  agree 
quite  closely  with  the  exact  (numerical)  pole  positions  shown  in  Fig.  9.  However,  while  the  numerical 
pole  calculation  only  furnishes  the  n -value  of  each  pole,  use  of  (20b)  in  addition  determines  its  /-value, 
i.e.,  it  tells  us  which  poles  belong  to  a  layer  of  given  /  and  thus,  give  rise  to  the  /th  creeping  wave.  The 
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Fig.  9.  Poles  of  TM  modes  of  oscillations  of  a  conducting  sphere  in  the  complex  frequency  plane  loi  =  is).  Dashed  curves  connect 
zeros  of  given  n ;  solid  curves  connect  zeros  of  given  I  whose  residues  synthesize  the  fth  creeping  wave  (from  T9]). 

asymptotic  expression  of  (20b)  thus  serves  to  uniquely  and  rigorously  establish  the  correct  connection 
between  SEM  poles  and  creeping  waves. 

Note  that  the  (solid-curve)  pole  layers  in  Fig.  9  that  belong  to  the  /th  creeping  wave  are  non -symmetric 
about  the  horizontal  axis.  The  statement  of  [7],  namely  that  for  reasons  of  the  overall  reality  of  the  field 
solutions  to  the  scattering  problem,  the  poles  have  to  be  symmetrical  about  this  axis,  is  not  contradicted 
thereby  since  inspection  of  Fig.  9  shows  that  the  individual  poles  are  located  symmetrically  above  and 
below  the  horizontal.  Only  their  ordering  into  creeping-wave  layers  is  non-symmetric.  Evidently,  individual 
creeping  waves  may  be  described  by  complex  fields,  but  the  total  field  obtained  by  summing  over  all 
creeping  waves  is  real.  Note  also  that  for  the  case  of  an  infinite  cylinder,  the  individual  SEM  poles  were 
found  [36]  to  be  symmetrical  again,  but  here  the  situation  is  more  complicated  since  has  a 

branch  cut  running  from  the  origin  of  the  s -plane  horizontally  to  the  left,  dividing  the  upper  and  lower 
s-plane  into  two  different  sheets  of  which  only  those  sheets  containing  symmetrical  poles  appear  to  be 
the  physical  ones.  No  branch  cut,  and  thus  only  one  single  sheet,  exists  for  the  present  case  of  the  sphere. 

The  SEM  poles  can  be  exhibited  in  the  denominator  of  (19)  by  Taylor-expanding  [jtfiL"(x )]’  about  its 
zeros.  Calling 

Ti(6)  =  (1/fcr)  ei,<tr_“"'sTM  (6), 


(21a) 


K~ 
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one  obtains 
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i ikr-uit  \ 


Ti(0)  =  - 


\r/a 


X  X 

1  =  l  n  -  1 


(-1  )'■ 

n+i  [xnij„(xn,)]' 

1  2 

dP„(cos  6) 

x  -X„, 

Xnl  [x„ih'„u Li„i)]" 

\n(n  + 1) 

d6 

(21b) 


for  backscattering  (0  =  0),  the  last  factor  becomes  unity.  The  first  factor  after  the  summation  sign  represents 
the  SEM  poles.  Proceeding  now  to  the  transient  case,  we  assume  an  incident  5-  pulse,  5(2  -4  ct ),  and  with 
a  dimensionless  time  variable  t  =  ( ct-r)/a  we  obtain 


Tx(6,  r)=  X  TT(0,t) 
i=i 


where  for  8  =  0: 


7T(0,r)=  —  I  (-1) 

r  n  =  1 


„  o  + 1  [Xnlin(Xnl)]'  -Unlr 

Xnl 


(21c) 


(2 1  d) 


In  view  of  (8),  this  expression  has  the  standard  SEM  form  for  transient  radar  echoes,  i.e.,  it  consists  of 
a  superposition  of  damped  sinusoids  (sometimes  referred  to  as  a  Prony  series).  While,  however,  SEM 
theory  applied  this  expression  to  (21c)  for  T\{0,  r)  as  a  whole,  we  consider  in  (21d)  separately  the  terms 
T{'\0,  t),  i.e.,  a  Prony  sub- series  summed  only  over  poles  in  one  layer  of  given  l,  the  reason  being  that 
this  subseries  synthesizes  the  /th  creeping  wave. 

In  Fig.  10  we  plot  on  an  arbitrary  scale  the  quantity  A'(r)  =  rRe  r'i"(0,  r),  i.e.,  the  backscattering 
amplitude  for  the  first  creeping  wave.  It  has  here  been  evaluated  [37]  using  the  asymptotic  expressions 
for  the  poles  of  (20b).  While  the  individual  damped  sinusoids  in  (2  Id)  are  quite  large  numerically  especially 
near  r  =  0,  they  all  cancel  in  the  sum  except  near  t  =  tt,  r,  etc.,  which  are  the  correct  arrival  times  of 
a  multiply-circumnavigating  creeping  wave  with  velocity  c.  The  exact  arrival  times,  however,  are  deter¬ 
mined  by  the  group  velocities  of  the  creeping  waves.  Using  the  fact  that  asymptotically,  P\( cosfl)oc 
cos[(n  +2)8  ~(tt/4)],  and  inserting  this  and  (20b)  in  (21d),  one  may  evaluate  the  ensuing  sum  by  the 
method  of  stationary  phase.  It  can  be  shown  that  the  ratio  of  spherical  cylinder  functions  asymptotically 
has  no  exponential  behavior.  If  to  lowest  order,  x„t  is  approximated  by  n  +j,  the  sum  has  stationary 
phase  for  r  =  (2N  +  l)ir  where  N  =0, 1,  2, ... ,  proving  the  multiple  (TV- fold)  circumnavigations.  To  higher 
order,  the  phase  <b„t  of  (2 Id),  i.e., 


rf)„i  =  (2TV+  1)-jt«  -xnrr 


(22) 


(with  x„t  =  Re  xnt),  must  be  stationary,  i.e.,  its  derivative  with  respect  to  n  must  vanish.  Noting  that  phase 
and  group  velocities  of  the  creeping  waves  are  given  by  [8] 


Ci(x)/c  =x/Re[ri(U)  +  !], 
crU)/c  =  l/Re[dn((jc)/djc], 


(23a) 

(23b) 


one  recognizes  using  (20a)  that  the  stationary-phase  condition  on  (22)  leads  to  the  dimensionless  arrival 
time  of  the  backscattered  pulses: 


rarr  =  (2N+lW(c?7c), 
or,  going  back  to  actual  time: 

t„„  =  (2/V  +  1  )-rra/cfp, 


(23c) 


(23d) 
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Fig.  10.  Pulses  of  /  =  1  (TM)  creeping  wave  in  the  scattering  amplitude  -r  Re  T\  (0,  r)  as  synthesized  by  the  appropriate  Prony 
subseries,  evaluated  by  the  stationary-phase  method  (from  [37]),  and  plotted  vs.  a  new  time  variable  a  *  r/[(2N  +  1  )ir  +  6],  as 

measured  from  the  individual  pulse  arrival  times  at  a  =  1  onward. 

showing  that  the  creeping-wave  pulses  propagate  with  the  group  velocity.  The  actual  pulse  shapes  of  the 
first  TM  creeping  wave,  shown  in  Fig.  10,  are  those  obtained  at  or  shortly  after  (*s5°)  one-half  ( N  =  0, 
top),  1^  (JV-1,  center)  and  2\  (N  =  2,  bottom)  circumnavigations  of  the  sphere,  by  replacing  the  sum 
over  n  in  (2  Id)  by  an  integral  which  was  evaluated  using  the  stationary  phase  method.  The  flowing-apart 
of  the  pulse  due  to  dispersion,  and  its  attenuation  in  the  course  of  the  multiple  circumnavigations,  due 
to  radiation,  is  evident. 


6.  Analysis  of  transient  creeping  wave  experiments 

Experiments  on  the  scattering  of  short  electromagnetic  pulses  from  conducting  spheres  and  cylinders 
have  been  carried  out  by  B.Z.  Hollmann  of  the  Naval  Surface  Weapons  Laboratory,  Dahlgren,  VA, 
USA  [38].  The  experimental  arrangement  consisted  of  a  conducting  ground  plane  with  a  vertical  conical 
transmitter  (S)  in  the  center  whose  tip  touched  the  ground  plane;  its  length  was  13  ft  9  in.  and  its  half 
opening  angle  5.29  degrees,  and  it  generated  a  near-plane  wave  field  with  a  vertical  electric  vector.  The 
receiver  antenna  (R),  a  half  TEM  horn,  rested  on  the  ground  plane,  as  did  the  targets:  a  hemisphere  as 
well  as  a  lying-down  half-cylinder  whose  mirror  image  at  the  plane  completed  them  to  a  full  sphere  and 
cylinder,  and  an  upright  cylinder  whose  length  was  doubled  by  the  mirror  image. 
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The  transmitter  generated  a  1-ns  pulse  with  a  smaller  negative  overshoot  of  1.5  ns  duration.  It  could 
be  represented  by 


P.ncU)=  Gftolexpf-iw/idw, 

J  -00 


with  a  spectrum  G(cj)  determined  from  the  experimental  pulse  shape  by  Fourier  inversion.  We  shall 
here  discuss  the  experimental  case  of  the  half-cylinder  lying  flat  on  the  ground  plane,  shown  in  top  view 
in  Fig.  11.  The  echo  pulses  then  consisted,  in  their  order  of  arrival  time,  of  a  specularly  reflected  pulse 
(1),  two  coinciding  pulses  (2,2')  reflected  from  the  cylinder  edge,  and  a  helical  creeping-wave  pulse  (3). 


Fig.  11.  Geometry  ot  radar  pulse  scattering  experiment  observing  helical  creeping  waves  (3)  on  a  conducting  cylinder. 


The  reflected  pulse  shapes  are  close  replicas  of  the  incident  pulse,  with  partially  calculable  distortions 
[24].  The  creeping  wave  pulse,  however,  is  distorted  due  to  the  dispersion  contained  in  (16c).  Its  amplitude, 
after  having  traveled  a  distance  s  =  (a202  +  z2)1/2  over  the  cylinder  surface,  is  given  by 


Percept,  0  =  |  GMexp  i«^  -  tj  dw 


where  to  lowest  order,  o  may  be  taken  from  (16c)  since  TE  creeping  waves  on  conducting  cylinders  and 
acoustic  creeping  waves  on  rigid  cylinders  have  the  same  phase  velocity,  to  this  order  [39]  (including 
the  value  of  q i).  It  turned  out,  however,  that  due  to  the  width  of  G(«),  a  more  accurate  representation 
was  needed.  Accordingly,  a  five-term  representation  of  nt  for  the  cylinder,  analogous  to  (16b),  was  used 
as  given  by  Franz  and  Galle  [39]. 

Figure  12  shows  the  measured  signal  vs.  time  corresponding  to  the  experiment  of  Fig.  11,  after 
subtraction  of 

(i)  instrumental  clutter  (measured  without  target), 

(ii)  the  specularly  reflected  pulse,  and 

(iii)  the  edge-diffracted  pulses. 

While  the  distortion  of  (ii)  could  be  taken  into  account  accurately  [24],  (iii)  is  a  hard  to  calculate  for  a 
curved  edge  and  we  subtracted  a  replica  of  the  incident  pulse  instead,  resulting  in  an  incomplete  subtraction 
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Fig.  12.  Experimental  helical  creeping-wave  pulse  on  a  conducting  cylinder,  compared  with  calculated  dispersive  (solid  curvel  and 

non-dispersive  (broken  curvel  creeping  wave  pulses. 

in  the  region  of  1-5  ns  after  the  specular  arrival  (the  original  height  of  the  latter  amplitude  was  ~80  on 
the  scale  of  Fig.  12). 

Overlaid  on  the  figure  is  the  creeping  wave  calculated  from  (24b)  (solid  curve),  as  well  as  an  artificial 
creeping  wave  calculated  assuming  no  dispersion  (dashed  curve),  which  is  thus  an  exact  replica  of  the 
incident  pulse.  The  agreement  with  the  observed  creeping  pulse  is  striking,  in  view  of  the  fact  that  no 
free  parameters  entered  the  calculation.  The  agreement  thus  includes  the  arrival  time,  polarity,  amplitude, 
and  dispersion  of  the  creeping  waves.  Similar  calculations  were  carried  out  for  the  sphere,  and  also  for 
upright  cylinders  (TM  modes).  Due  to  the  higher  attenuation  of  these  modes,  no  creeping  waves  could 
be  observed  in  the  latter  case,  however. 


7.  Resonance  scattering  from  targets  of  arbitrary  shape:  Acoustic  case 

While  the  above-discussed  examples  treated  analytically  the  scattering  resonances  in  targets  of  separable 
geometry,  bodies  of  arbitrary  (non-separable)  shape  are  obviously  of  very  great  importance.  Field 
expansions  in  spherical  coordinates,  such  as  (la,  b),  (12)  or  (14a),  can  still  be  used  for  this,  but  the 
non-separable  shape  of  the  target  will  induce  a  coupling  of  the  channels,  i.e.  the  nth  normal  mode  of 
the  incident  wave  will  be  coupled  to  all  normal  modes  of  the  scattered  wave,  and  vice  versa.  The  T- 
function  Tn  of  (14),  e.g.,  will  thus  be  replaced  by  a  ‘T  matrix'.  The  corresponding  T  matrix  method  was 
introduced  by  Waterman  [40]  (see  also  Pao  [41]),  and  we  shall  use  it  as  a  basis  of  our  further  discussion. 

The  T-matrix  method  has  been  combined  by  us  with  the  resonance  approach  in  the  acoustic  case  [42]. 
Originally,  Waterman  had  also  applied  his  T-matrix  approach  to  radar  scattering  from  conducting  bodies 
[43],  and  we  have  adapted  his  computer  program,  which  originally  was  designed  for  the  calculation  of 
scattering  amplitudes  only,  to  the  task  of  obtaining  the  complex-frequency  (SEM)  poles  of  the  scattering 
amplitudes  for  conductors  of  arbitrary  shape  [44], 

Examples  of  spheroids  and  finite-length  cylinders  have  been  treated  in  this  way,  and  we  have  furthermore 
generalized  Waterman’s  scattering  code  which  was  originally  restricted  to  the  case  of  axially-symmetric 
bodies  [43],  to  targets  of  nonaxial  shape  so  that  the  SEM  poles  of  ellipsoidal  targets  could  be  obtained 
also  [44]. 


H.  Uberall  et  al.  /  Resonance  scattering  theory 


2  1  prolate 

SPHEROID 


Fig.  13.  Shifts  and  splittings  of  SEM  poles  in  the  complex  frequency  plane  upon  deformation  of  a  conducting  sphere  into  a  2: 1 

prolate  spheroid  (from  [23,44]). 

The  following  two  figures  present  some  of  our  results  on  these  topics.  Figure  13  plots  the  (external) 
TM  poles  of  a  conducting  sphere,  identical  to  those  in  Fig.  9,  as  large  circles  labeled  by  (nMl),  as  well 
as  TE  poles  labeled  by  (nEl),  and  shows  how  they  shift  and  split  as  the  sphere  is  deformed  into  a  2 : 1 
spheroid  by  stretching  the  axis  of  symmetry  (b),  while  the  short  axes  are  kept  equal  to  a.  The  split 
components  are  labeled  by  an  azimuthal  quantum  number  m.  Previous  work  [7]  had  only  obtained  the 
m  =  0  components  and  their  shifts,  but  had  ignored  the  splitting  of  the  poles  into  various  m  components. 

A  different  way  of  presenting  these  results  can  be  taken  by  plotting  the  real  parts  of  the  pole  frequencies 
only,  in  the  form  of  a  ‘level  digram’  familiar  from  atomic  physics.  This  is  done  in  Fig.  14,  and  it  seems 
proper  calling  this  approach  'Radar  Spectroscopy'  [23],  in  analogy  to  optical  spectroscopy  which  is  based 
on  atomic  level  diagrams.  Figure  14  shows  the  ( nEl)m  spectra,  for  /  =  1  (the  least-damped)  creeping-wave 
poles  only,  successively  for  a  conducting  sphere,  a  1 : 1  cylinder,  a  2  : 1  spheroid,  a  2 : 1  cylinder  and  a 
very  long  (100: 1)  cylinder,  the  latter  taken  from  [7].  The  point  of  such  a  radar  spectral  diagram  is  that 
it  indicates  the  sensitivity  of  the  levels  to  changes  in  shape  of  target,  hence  offering  the  possibility  of 
target  identification  and  the  determination  of  target  shape  from  an  experimentally  obtained  resonance 
spectrum  of  the  target.  This  point  will  be  elaborated  on  below. 


8.  Acoostic  spectroscopy:  The  inverse  problem 

An  ‘acoustical  spectrogram'  similar  to  that  of  Fig.  14  has  been  obtained  by  us  for  the  acoustical  case 
[22],  and  is  shown  in  Fig.  15.  It  presents  the  (internal)  real,  normalized  eigenfrequencies  x  =u>a/c  of 
fluid  objects  in  vacuum  (which  is  a  good  approximation  for,  e.g.,  water  drops  in  air),  and  shows  the  level 
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Fig.  14.  Radar  Spectroscopy'  diagram  of  the  real  parts  of  SEM  pole  positions  of  a  conducting  sphere,  1  : 1  cylinder,  2  : 1  spheroid, 
2 : 1  cylinder,  and  l(M):  1  cylinder  </  =  1  TE  poles  only),  showing  level  shifts  and  splittings  (from  [23]). 


shifts  and  splittings  as  the  sphere  is  deformed  into  a  b/a  =  1 . 1 1  . . . :  1  prolate  spheroid  or  circular  cylinder, 
a  4: 1  spheroid  and  cylinder,  and  finally  an  infinite  cylinder.  The  levels  are  degenerate  with  respect  to 
azimuthal  quantum  number  m  for  the  sphere,  then  split  up,  but  become  degenerate  again  with  respect 
to  the  axial  quantum  number,  for  the  infinite  cylinder. 

These  eigenfrequencies  were  here  determined  by  us  not  via  a  use  of  the  T-matrix  method,  but  by 
searching  for  the  zeros  in  frequency  of  radial  spheroid  eigenfunctions  on  the  boundary;  the  internal 
cylinder  problem,  however,  is  soluble  in  an  elementary  fashion. 

If  experimentally,  a  sufficient  portion  of  the  eigenfrequency  spectrum  of  the  target  has  been  determined 
e.g.  from  an  analysis  of  the  scattered  echoes,  an  attempt  at  target  classification  (i.e„  solution  of  the 
inverse  problem),  can  be  made.  For  the  sphere,  the  eigenfrequencie*  are  determined  from  the  zeros  of 
the  spherical  Bessel  functions,  and  if  the  measured  spectral  level  ratios  agree  with  the  Bessel  function 
zeros,  the  spherical  nature  of  the  target  is  established.  The  absolute  value  of  any  eigenfrequency  then 
provides  us  with  the  ratio  a/c  where  c  is  the  sound  speed  inside  the  target.  If  the  spectrum  shows  the 
sequence  of  spheroidal  levels,  the  spheroidal  nature  of  the  target  is  established  and  the  ratios  a/c  and 
b/c  can  be  determined  by  absolute  level  measurements. 
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Fig.  15.  'Acoustical  Spectroscopy'  diagram  for  the  eigenfrequencies  of  fluid  objects  in  vacuum:  sphere,  prolate  spheroids,  finite 
and  infinite  cylinders,  showing  level  shifts  and  splittings  (from  [22]). 


An  experimental  determination  of  the  eigenfrequency  spectrum  of  an  elastic  cylinder  by  acoustic 
scattering  has  recently  been  accomplished  by  Maze,  Taconet  and  Ripoche  [45]  in  an  ingenious  fashion. 
If  a  conventional  tone-burst  (near-steady  state)  experiment  is  performed,  the  amplitude  measured  at 
various  frequencies  consist  of  an  interfering  superposition  of  resonances  and  background  in  the  fashion 
of  (7d),  from  which  the  resonances  are  often  hard  to  separate,  due  to  the  fact  that  the  measured  data 
show  all  resonances  contained  in  the  sum  over  modes  n.  Part  (b)  of  Fig.  16  presents  such  a  measured 
spectrum  consisting  of  resonances  and  background,  for  the  case  of  an  aluminum  cylinder  in  water.  If, 
however,  one  waits  until  the  specularly  reflected  background  has  passed,  the  slightly  attenuated  Whisper¬ 
ing-Gallery  waves  which  keep  circulating  around  the  cylinder  in  a  resonant  fashion,  will  then  provide 
an  echo  signal  that,  upon  a  frequency  sweep,  shows  their  resonance  spectrum  (part  (a)  of  Fig.  16).  This 
procedure  represents  a  direct  experimental  way  of  determining  the  pure  acoustic  spectral  diagram  of  a 
target,  without  the  interference  of  a  refiected-wave  or  external  resonant  background,  hence  offering  a 
direct  approach  to  acoustic  spectroscopy  and  target  classification.  Indeed,  part  (c)  of  Fig.  16  shows  the 
calculated  (top)  and  measured  (bottom)  resonance  frequencies  of  the  cylinder,  labeled  by  n  (mode 
number)  and  /  (surface  wave  label).  Such  a  comparison  offers  the  possibility  of  a  determination  of  the 


H.  Uberali  el  at.  /  Resonance  scattering  theory 


327 


>. 


l 

15 

— 1 
20 

- T" 

25 

H - 

30 

1 

35  ko 

1 

b 

u 

b 

b 

|7 

18 

|9 

||0  III  2 

ii 

.12 

13 

u 

15 

16 

i? 

le  3 

h 

lz 

13 

.  ,4 

15 

.  16  4 

Ii 

12 

Is 

U  5 

h 

|2 

13  U  6 

Ii 

12  7 

ll  8 

”1 

15 

1 

20 

1 

25 

1 

30 

1  * 

35  ka 

b 

|4 

Is 

|6 

2 

h 

Iz 

13 

u 

|6 

I? 

3 

2 

13 

14 

b 

<6  4 

12 

5 

Ii 

6 

12  7 

20 


25 


30 


I 

35  ko 


Fig.  16.  Measured  acoustic  spectrum  of  an  aluminum  cylinder  in  water,  (a)  pure  resonances,  (b)  resonances  and  background,  (c) 
comparison  of  measured  and  calculated  resonance  frequencies  (from  [45]). 

material  composition  of  the  target  via  acoustic  spectroscopy,  just  as  optical  spectroscopy  determines  the 
chemical  composition  of  a  substance.  The  term  'Acoustic  Spectroscopy’  has  been  coined  in  this  connection 
by  Derem  [46], 


9.  Summary 

The  preceding  examples  illustrate  the  accomplishments  of  the  Resonance  Scattering  Theory  [1]  in  a 
great  variety  of  fields  of  classical  physics  (radar,  sonar,  elastic  waves  and  geophysics)  regarding  two 
principal  aspects,  namely: 

(i)  providing  analytical  connections  between  the  resonances  of  a  target  (the  ‘SEM  poles’)  and  the 
(internal  and  external)  surface  waves  on  the  target  (described  by  'Regge  poles’),  and  hence  establishing 
a  physical  interpretation  of  the  resonances,  and 

(ii)  providing  a  potential  means  for  the  solution  of  the  inverse  scattering  problem,  i.e.  the  identification, 
and  the  material  and  shape  determination  of  the  target  from  an  analysis  of  the  resonance  effects  contained 
in  the  returned  echo  signals. 

In  relation  to  the  inverse  problem,  we  have  here  proposed  to  represent  the  spectrum  of  resonance 
frequencies  of  a  given  target  by  a  level  scheme'  similar  to  that  of  atomic  spectra,  and  have  accordingly 
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termed  this  approach  'Radar  Spectroscopy'  or  ‘ Acoustic  Spectroscopy'  (see  also  [46]).  The  shifting  and 
splitting  of  the  levels  during  the  changes  of  shape  of  a  target  of  known  composition  (and,  by  implication, 
the  analogous  level  shifts  undergone  during  changes  of  composition)  which  we  have  illustrated  here,  then 
offer  a  possibility  for  target  identification,  and  hence  for  a  solution  of  the  inverse  problem.  Some  experience 
will  be  required  when  putting  such  a  scheme  into  working  order,  as  it  has  been  the  case  in  optical 
spectroscopy  as  well.  In  addition  to  merely  the  level  positions  and  their  widths,  the  intensities  of  the 
resonance  excitations  (i.e.,  the  SEM  pole  residues)  will  provide  crucial  information  on  these  and  on  other 
aspects  of  the  inverse  problem  such  as  target  orientation,  as  will  the  use  of  polarization  data.  (Note  that 
level  positions  and  resonance  widths  are  inherent  properties  of  a  target  while  their  excitation  strengths 
depend  on  the  means  of  excitation  chosen,  such  as  the  incident  beam  direction  and  polarization  with 
respect  to  the  target  geometry.) 

From  the  foregoing  discussion,  there  emerges  a  global  picture  of  target  resonance  effects  that  promises 
to  result  in  a  useful  tool  for  target  identification  and  for  a  solution  of  the  inverse  problem.  The  mentioned 
recent  experimental  work  [45],  which  has  provided  a  way  of  exhibiting  the  pure  scattering  resonances 
in  a  direct  fashion  for  the  acoustic  case,  should  be  particularly  relevant  in  this  respect. 
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